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i—H ■ In many inflationary models, a large amount of energy is transferred rapidly 

to the long-wavelength matter fields during a period of preheating after in- 
flation. We study how this changes the dynamics of the electroweak phase 
transition if inflation ends at the electroweak scale. We simulate a classi- 
cal SU(2) xU(l)+Higgs model with initial conditions in which the energy is 
■ concentrated in the long-wavelength Higgs modes. With a suitable initial en- 

Q^ ■ ergy density, the electroweak symmetry is restored non-thermally but broken 

again when the fields thermalize. During this symmetry restoration, baryon 

t— ( \ number is violated, and we measure its time evolution, pointing out that it is 

highly non-Brownian. This makes it difficult to estimate the generated baryon 

^5 asymmetry. 
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I. INTRODUCTION 



In order to explain the observed baryon asymmetry of the universe, a theory must sat- 
isfy three conditions [0; baryon number violation, C and CP violation and deviation from 
thermal equilibrium. Grand unified theories offer the required properties, but in addition 
to baryons, they also lead to formation of magnetic monopoles. This monopole problem 
can be solved by a subsequent period of inflation, but it would also wipe out any baryon 
asymmetry. 

Although the baryon number is conserved perturbatively in the electroweak theory, there 
are non-perturbative effects that violate it. At zero temperature, this effect is negligible, 
because the baryon number can only be changed by strongly suppressed tunneling processes 
0. However, there are classical field configurations, sphalerons M, that interpolate between 
different baryon numbers and become unsuppressed at high temperatures . Thus, provided 
that CP is violated, the electroweak theory can explain the baryon asymmetry by itself, as 
long as the electroweak phase transition gives the necessary non-equilibrium state. 

However, sphaleron processes do not stop immediately after the electroweak transition, 
and thus the baryon asymmetry can be washed out by the sphalerons, unless the electroweak 
symmetry is already strongly broken when the universe returns to thermal equilibrium. In 
the standard big bang scenario, only a strongly first-order phase transition can prevent 
this baryon washout, but lattice Monte Carlo simulations have shown that the thermal 
transition cannot be strong enough in the standard model |5[]. Although it is possible in 
some extensions, such as MSSM ||, this has led people to consider other possibilities. 

A particularly attractive idea, suggested in Refs. 0,H, is based on the non-equilibrium 
dynamics at the end of inflation. When inflation ends the universe is very cold and all the 
energy is stored in the inflaton field, which subsequently decays into the standard model 
fields, thus heating up the universe to the reheat temperature T rh . Instead of ordinary 
perturbative decays this can take place non-perturbatively, and much faster, if the inflaton 
starts to resonate with the matter fields ||. In this period of parametric resonance, called 
preheating, the energy is transferred mainly to the modes with longest wavelengths, and 
they reach an effective equilibrium state at a temperature T e g which can be much higher 
than T r h- This makes it possible to effectively restore symmetries that are spontaneously 
broken at T r h ]TI]|JTT[ |. Eventually, the universe thermalizes and cools down to the equilibrium 
temperature T r h, undergoing a non-thermal phase transition to the broken phase. 

In the special case in which T r h is below the electroweak critical temperature, baryon 
number is violated by sphalerons and the system is far from equilibrium during the period 
of non-thermal symmetry restoration, and therefore the conditions for baryogenesis are sat- 
isfied MM . The washout problem is avoided if T rh is low enough, because the rate at which 
the universe cools is determined by the decay rate of gauge bosons, and is much faster than 
the decay rate of baryons. 

Because of the non-perturbative, non-equilibrium nature of the process, the only reliable 
way of studying it is to simulate it numerically. Also, the fluctuations produced during 
preheating have large occupation numbers so they can be considered as interacting classical 
waves |T2]JT3|] . Thus the dynamics of the system can be studied simply by solving the classical 
equations of motion numerically, which is a relatively straightforward task. 

Non-thermal symmetry restoration was first observed in numerical simulations in a scalar 
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field model in Ref. More recently, two of us studied the same process in the Abelian 
Higgs model j0| , which has a very similar structure to the electroweak theory. The results of 
the simulations support the scenario, but unfortunately there is no analogue to the Chern- 
Simons number in the three-dimensional Abelian Higgs model, and thus the question of 
baryogenesis could not be addressed directly. The generated baryon asymmetry has been 
measured in one-dimensional Abelian simulations ||, but it is not possible to convert the 
results into three dimensions, and therefore a quantitative test of the scenario is still missing. 

In this paper, we study the dynamics of the full electroweak theory in non-equilibrium 
conditions. The methods we develop can be applied to studies of baryogenesis during pre- 
heating, in a given inflationary model. However, such simulations would necessarily be very 
model-dependent and therefore we consider only a simplified scenario in this paper. In 
that scenario, we find that the rate of true topological sphaleron transitions is slower than 
expected because of the interaction between the gauge and Higgs fields. 

We start by reviewing the scenario of electroweak baryogenesis at preheating in Sect. [TI]. 
We only simulate the CP invariant case, but since CP violation is necessary for baryogenesis, 
we discuss the possibilities of including it in the simulations in Sect. In Sect. |Vj, we 
discuss the approximations that are needed in order to make the simulations feasible and 
show that all of them are justified. Numerical simulations also require that the space-time 
is discretized, and we present the details of this in Sect. [V|. In Sect. [VI|, we discuss how 
the requirement that the baryon number is not washed out when the system thermalizes, 
restricts the parameters of the theory. In Sects. [VI I| and |VIII| we present the results of our 
simulations, and in Sect. |X|, we discuss the conclusions we can draw from them. Some 
technical details about the simulations are presented in the appendices. 



II. BARYOGENESIS AT PREHEATING 

The scenario of electroweak baryogenesis at preheating |7]J§ assumes a model of inflation 
in which the energy scale of the inflation is so low that the equilibrium temperature of the 
universe never exceeds the electroweak critical temperature after inflation. The inflaton is 
assumed to couple to the standard model Higgs field. 

The inflation dilutes the inhomogeneities in all the fields away and results in a cold, 
vacuum configuration. However, the zero mode of the inflaton has a large value, and after 
inflation ends it oscillates with some frequency to. This starts a parametric resonance with 
those Fourier modes of the Higgs field that have the same frequency [0]. This resonance 
amplifies the quantum fluctuations of the Higgs field exponentially and transfers very rapidly 
a significant fraction of the energy of the inflaton to the Higgs field. If the frequency u is small 
enough, the resonance only excites the long-wavelength part of the Higgs power spectrum, 
but the occupation numbers of those modes becomes huge. 

The energy subsequently flows into the gauge bosons because of their strong coupling to 
the Higgs field, and soon a quasi-equilibrium state is reached in which the bosonic modes 
below a certain cutoff are approximately in thermal equilibrium while the other modes are 



still in vacuum |13|. Since the energy density is concentrated in a small number of degrees of 



freedom, the effective temperature T e g of the long-wavelength modes is much higher than the 



reheat temperature T r h, and therefore the electroweak symmetry can be restored [10,[IT,15 
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We can estimate that k* ~ awT e s |15j, because above that scale the interactions are per- 
turbative and therefore the thermalization rate of the modes with k > a.wT e & is suppressed 
by powers of coupling constants. It is important that the scale k* is the same as the typical 
size of a sphaleron configuration ||. Thus the sphaleron transitions are insensitive to the 
physics at higher momenta, and we can try and estimate their rate in the same way as in 
thermal equilibrium. In the broken phase, the sphaleron rate r sp h, defined as the number of 
sphaleron processes per unit volume per unit time, is exponentially suppressed by the Higgs 
expectation value, but in the symmetric phase it is unsuppressed. In the presence of an 
ultraviolet cutoff at temperature T eS , the classical sphaleron rate r sph in the symmetric 



phase has been estimated to be ]I6| 



sph » a 5 w T^/K ~ c4T c 4 ff . (1) 
Each sphaleron process changes the SU(2) Chern-Simons number 

iVcs(i) - iVcs(O) = J* dt J d 3 xe^W* u W; a (2) 
by one, leaving the U(l) equivalent, 

ncs(t) - ^cs(O) = fdt [ d 3 xe^B^B p(T (3) 

OZ7T z Jo J I 

unchanged. Therefore, Nqs changes wildly during the non-thermal symmetry restoration, 
and the baryon current is not conserved, due to the anomaly, 



327T i 



(4) 



The system is strongly out of equilibrium and therefore baryon asymmetry can be generated 
provided that CP violation is strong enough. As was pointed out in Ref. [17|], the oscillations 
of the inflaton can enhance this process. 

The effective temperature T e g decreases gradually, partly because the bosonic modes with 
higher momenta become excited but mainly because of decays of the bosons into fermions. 
The experimental value for the decay rate, 7^2 GeV, is much higher than the rate of 
baryon number violation 

IB - f ^, (5) 

which implies that the sphaleron processes stop instantaneously when the system enters the 
Higgs phase, and the baryon asymmetry generated during the period of symmetry restoration 
is not washed out. Simulations in a 1+1D toy model support this scenario ||, but obviously 
they cannot decide whether it can really work in the context of the realistic standard model. 
Our aim in this paper is to partially address that question. 

We are not going to address the issue of finding a model of inflation that would have a 
low enough reheat temperature and would produce the necessary initial conditions. In fact, 
we don't treat the inflaton as a dynamical field, but simply assume that the state of the 
system after the parametric resonance has transferred the energy to the Higgs field, can be 
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approximated by a configuration in which the Higgs has a large expectation value and its 
inhomogeneous modes are in vacuum. This means a drastic simplification in comparison 
with the scenario proposed in Refs. and is obviously only a first approximation, but it 
allows us to parameterize the properties of the inflaton by a single number (p Q , the Higgs 
expectation value after the preheating, and therefore makes the analysis simpler. In any 
given inflationary model, it is a straightforward task to include the inflaton as a dynamical 
field, thus making the treatment much more accurate, but at the expense of making the 
results specific to that model. 



III. CP VIOLATION 

In the absence of CP violation, there is no preferred direction for the change of the baryon 
number, and therefore domains of matter and antimatter would be formed. In practice, these 
domains are far too small to explain the baryon asymmetry of the universe, and therefore 
CP violation is necessary. Unfortunately, the only experimentally confirmed source of CP 
violation in the standard model, due to Cabibbo-Kobayashi-Maskawa mixing of quarks, is 
too small to cause a sufficient baryon asymmetry |18|JT9|] . This means that we must resort 
to extentions of the standard model to find additional sources of C and CP violation. 

Typically, effects of CP violation are modelled through an effective field theory approach. 
Once all degrees of freedom except the gauge fields, the Higgs, and the inflaton have been 
integrated out, the effective Lagrangian will contain non-renormalizable operators that break 
CP. The lowest of these is a dimension-six operator JT9I 



where the dimensionless parameter 5 cp is an effective measure of CP violation, and M new 
represents the scale at which the new physics, responsible for this effective operator, is 
important. Assuming M new ~ 1 TeV, the experimental constraints on the magnitude of the 
neutron electric dipole moment |2(J lead to an upper bound 5qp < 1, for the CP violating 



parameter ||21|| . Recent experimental constraints emerging from measurements of the electric 
dipole moment associated with mercury, indicate that the bound on Sqp could be much 
tighter f2"2]] , being reduced to 5cp < 0.05. Such a strong bound has to be taken with a 
degree of caution as its derivation assumes a number of things about the magnitude of the 
contributing sources. 

There are two ways that we know of to calculate the baryon number when CP violation 
is present, each with their own problems. Firstly one could perform numerical simulations 
with an explicit CP violating term in the equations of motions, as was done in H. However, 
as explained in [p3fl , such terms also break the lattice gauge invariance which is a rather high 



price to pay. Although not an insurmountable problem [23] we shall not tackle this here. 
The second approach is to treat 5cp perturbatively and not evolve the system with the CP 
violating term, rather use this evolution as a background. For this one measures the change 
in the Chern-Simons number and relates this, via a Boltzmann equation, to the change in 
the baryon number. This approach also does not work here as it relies on Acs following a 
random walk which, as we shall argue later, is not the case here, so we do not attempt to 
measure the baryon asymmetry in this paper. 
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IV. APPROXIMATIONS 



Given the non-perturbative and non-equilibrium nature of the scenario, it is obviously 
not possible to study it in the full standard model without any approximation. Fortunately, 
it is possible to justify three approximations, which together make the problem tractable: 

• Bosonic approximation: Initially, most of the energy is in the Higgs field, and the 
parametric resonance will transfer it very effectively into the SU(2) and U(l) gauge 
fields. Because a similar parametric resonance cannot take place for the fermionic 
fields and because the Higgs couples relatively weakly to most of them, we assume 
that the fermions remain in vacuum. In that case, the main effect of the fermions 
is to give an extra channel into which the bosons can decay. We approximate this 
by adding a simple damping term with the damping rate 7 ~ 2 GeV, which agrees 
with the experimental decay with of W and Z bosons. This approximation does not 
take into account relativistic effects or the oscillations of the Higgs condensate, but 
we believe that it gives the correct qualitative picture. Improved accuracy could be 
obtained by introducing extra degrees of freedom that would mimic the effect of the 
fermions, but at the price of extra computational cost. 

Of course, these decays also transfer energy into the fermionic sector, which makes 
our approximation break down at t > 7 _1 . However, this does not restrict our analysis 
much, since if the transition to the Higgs phase has not taken place by this time, the 
reheat temperature will be close to the critical temperature, and therefore the baryon 
number is likely to be washed out. 

• SU(2)xU(l) approximation: The Higgs and the SU(2)xU(l) gauge fields interact 
with SU(3) gluons only indirectly via fermions. Therefore, if we assume that the 
fermions are in vacuum, the same will apply to the gluons as well. In fact, the flow 
of energy from fermions to gluons decreases the energy in the fermionic sector, and 
therefore makes the bosonic approximation remain valid for longer times. 

• Classical approximation: After preheating, the long-wavelength modes of the Higgs 
and the gauge fields have a large occupation number, and therefore they can be treated 



as classical fields ||12|jl3| . Consequently, time evolution of the system is given by the 
classical equations of motion, which can be solved numerically in a straightforward 
way. 

As a result of these approximations, the time evolution of the system is given in the 
temporal gauge by the classical equations of motion 

d 2 <f) = D t D4 + 2A [r, 2 - 0V) - 7<%0, 
d 2 Bi = -djBij + g'lm^Dicj) - ^d B u 



idoWu (7) 



where is the Higgs field, B t is the U(l) gauge field, Wi is the SU(2) gauge field, and the 
covariant derivative is 
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Di = di- ~gWi - l -g'B,. (8) 

Assuming that the decay rate 7 is the same for all the fields means that Gauss's law will 
be conserved. Thus it gives a useful check for the numerics of the simulation. If we define 
the canonical momenta 

7T = 9 O 0, E i = -d B il F^-doWi, (9) 
Gauss's law has the form 

d{Ei = g'lm^cf), 



[D h Fi] =ig 



7T0t — — <^7T — h.C. 



(10) 



The temporal gauge leaves the system invariant under time-independent gauge transfor- 
mations. This remaining gauge degree of freedom is fixed by choosing that the Higgs field 
is initially real and positive and that only its lower component has a non-zero value. 

Even though we use classical equations of motion, it is essential to include the quantum 
vacuum fluctuations in the initial configuration, because they act as seeds for the amplifi- 
cation of the low- momentum gauge field modes |12|]. They are approximated by classical 



fluctuations with the same equal-time two-point correlation as in the quantum theory at 
tree level. For each real field component Q of mass m and canonical momentum P, this 
means 

(Q*(t, k)Q(t, k')) = 1 (27r) 3 ^(fc - k'), 
2y k 2 + m 2 



(P%k)P(t,k')) = yk2 + m \ 2n)W\k-k>). (11) 
Since the Higgs is a complex field, it may be decomposed into four real components 

(12) 




Note that as <po has a large value so the masses of these components are different. 

We must also require that the initial configuration satisfies Gauss's law (|l0l) . This is 
achieved by a cooling procedure described in Appendix B. 



V. LATTICE 

In order to simulate the system numerically, we discretize it on a N 3 lattice with lattice 
spacing a in such a way that the Higgs field is defined at lattice sites and the gauge fields 
are represented by group elements at links between the sites. In the time direction, the time 
step is a t and the canonical momenta of the fields are defined at half-way between the time 
steps. The details of the discretization, as well as the lattice equations of motion, are given 
in Appendix A. 
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On a lattice, we label the lattice sites by integer labels Xj £ {!)•••> ^} an d the corre- 
sponding lattice momenta by Kj G {0, . . .N— 1}. We define the discrete Fourier transform 
by 

W) = ^3 E ex P f^Ett**/*) (13) 
The physical momentum vector is then given by 

ikj = - [1 - exp(-27iiKj/N)} (14) 

(X 

which has the property that the backward lattice derivative acting on a field in x space is 
equivalent to multiplying the k space field by ikj. The lattice Laplacian acting in x space 
also has the effect in k space of multiplying by —J2j \^j\ 2 , in analogy with the continuum 
case. 

The fields on the lattice are therefore required to satisfy [cf. Eq. (|TTD1 

1 1 

2^\k\ 2 + ™l N a 



Ejlkj^ + ml 1 



W*> = -^3%" ( 15 ) 

The lattice discretization regularizes the theory, and it must also be renormalized, i.e. the 
bare parameters used in the lattice equations of motion must be chosen in such a way that 
the values of any observables agree with those measured in experiments. To our accuracy, 
we can neglect the renormalization of the coupling constants, and therefore we only consider 
the Higgs mass counterterm. It can be determined by calculating the Higgs mass to one-loop 
order in lattice perturbation theory 

m% = mLt + (6A + \g 2 + \g' 2 ) (Q'Q) , (16) 

where the subscript in (. . .) indicates a free-field expectation value, and m 2 H is the physical 
Higgs mass. In the approximation where the masses inside the loops are taken to be zero, 
the correlation function (|i"5|) implies 



(Q*Q)o -E /Jv, / 1 * ^ A-2 r ^ f E si * 2 wj) 2 -^. (17) 

X ( Na ) 2JJ: Ak A 2 47r 3 a 2 7o W / a 2 



and m 2 H is the physical Higgs mass. This gives the bare mass parameter that must be used 
in the simulation, 

<tt « rn 2 H - [Q\ + -g 2 + -g' 2 j — — . (18) 

Also, we define the "renormalized" version of ($</>) by subtracting the one-loop diver- 
gence 
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non renormalized 
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FIG. 1. Evolution near the vacuum for a renormalized (bottom) and non renormalized (top) 
mass term with g = g' = 0.1 and a = 2,3,4 x 10~ 3 GeV^ 1 . The straight line is ry 2 . 



hatt — 2(Q*Q)o 



0.452 



/latt 



(19) 



This will be a useful quantity as an effective order parameter, but we emphasize that it 
is not a physical observable. In Fig. p] we start the fields near the vacuum solution and 
let them evolve. The upper graph shows how the equilibrium value of {(jr(j)) without the 
mass counterterm depends on the lattice spacing, relaxing to larger values as a is decreased. 
When the fluctuations are taken into account by the counterterm, lower graph, we see that 
this dependence is removed. Such behaviour is essential if we are to have results free from 
lattice artefacts. 

For our simulation we include a phenomenological damping term to account for the 
fermions. The effect of this is to damp the fluctuations on all length scales, this then 
changes the two point correlator ([15]) which appears in the mass counterterm flip]). When 
the counterterm is calculated one is interested in how the large momentum modes behave. 
In the limit k 3> 7 the dispersion relation becomes 



• 7 



(20) 



showing that (f)k(t) = exp(zu;fc£)0fc decays as e 7 *' 2 and thus the actual mass counterterm 
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cooling depth=0 

cooling depth=0.9a 2 

cooling depth=1 .8a 2 

cooling depth=4.5a 2 

— cooling depth=9.0a 2 




0.1 0.15 
time (GeV~ 1 ) 

FIG. 2. Measuring Nqs with different cooling paths. The path lengths are 0. 0.9a 2 , 1.8a 2 , 4.5a 2 , 
9.0a 2 . Without cooling, the ultraviolet noise blurs the signal, but with cooling the sphaleron can 
be seen at t = 0.1 GeV -1 . 



and (0V) counterterm (19) must be multiplied by e -7 *. 
We also need a lattice representative of the Chern-Simons number (|2]). There is no simple 
way to translate this continuum expression to a lattice, the essential problem being that any 



lattice definition of e^W'W* is not a total derivative [pcjfl . Here we follow the reasoning of 



Ambj0rn and Krasnitz |24j where a daughter set of fields are used to calculate eV^W^W". 
In fact, as pointed out in |25[ only the gauge fields should be included in this daughter set. 
This second set of fields are the cooled image of the first, removing UV fluctuations and 
leaving smooth gauge fields, bringing e^WIW* close to a total derivative. This method 
was developed in response to the realization that early numerical simulations of the Chern- 
Simons diffusion were measuring a large component of thermal noise rather than actual 
topological transitions. To test the method for measuring Acs we evolved the fields in a box 
small enough that only one sphaleron could fit in it, thus we should see only one sphaleron 
at a time and A C s should change by one unit. In Fig. |2| we display the difference such a 
cooling path can make in determining where a sphaleron transition has taken place. If we 
measure Acs without using the cooling path (upper curve) then the UV fluctuations mask 
the transition, but as the cooling depth is increased the integer transition between baryon 
vacua becomes clearer. In the following we shall use a cooling depth of 0.9a 2 . 



VI. PARAMETERS 

In the simplified model we are discussing here, the dynamics of the theory depends only 
on two parameters; the Higgs mass m# and the initial value 0o of the Higgs field. We 
choose the value m# = 100 GeV, but unlike standard electroweak baryogenesis, the scenario 
discussed here is probably not very sensitive to the Higgs mass. The coupling constants are 
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taken from electroweak data |26|]; g = 0.64, g' = 0.34, 17 = 250.7/^ GeV, 7 = 2 GeV. 

The initial value 0o of the Higgs field parameterizes the initial energy density, and there- 
fore also determines the final reheat temperature T r h to which the system eventually equi- 
librates. Although in the bosonic theory all the vacua with different Acs are degenerate, 
the one with zero baryon number has the lowest energy once the fermions are included 
and therefore the sphaleron transitions tend to wash out the baryon number. Thus, if T r h 
is above or too close to T c , the sphaleron rate at the final state can be so high that the 
baryon asymmetry generated before disappears. It has been estimated that to avoid this, 
one needs [[27j .|i"§ll 

v{T Th )>T A , (21) 
where v (T rh ) is the expectation value of the Higgs field. In one-loop perturbation theory, 



28| 



^P,= (l-gV. 1? = * (200 GeV)' (22) 

and thus the condition (^) becomes T r h < 155 GeV. The corresponding energy density is 
e(T r h) ~ (vT 2 /30)(7*T r j 1 , where « 100 is the effective number of degrees of freedom at the 
electroweak scale. Since e(T rh ) must be equal to the initial energy density e = A0q> we find 



30A \ 1/4 



rh 



O.220 o . (23) 



Consequently the constraint ( 21"D becomes 

0o < 700 GeV. (24) 



VII. SYMMETRY RESTORATION 

In Fig. ^ we show the time evolution of the "renormalized" value of (|0| 2 ) [see Eq. (fL9|)1. 
We use a 60 3 lattice with the initial condition O = 700 GeV [cf. Eq. and lattice spacings 
a = 3 x 10 -3 GeV -1 , a t = 0.2a; unless otherwise stated these will be used throughout the 
paper. We also performed test runs with smaller lattice spacings and time steps, finding 
results that were in statistical agreement with those presented here. 

Since (|0| 2 ) gets a positive definite contribution from every Fourier mode, it is not zero 
even in the symmetric phase. Furthermore, since there is no phase transition in the elec- 
troweak theory even in equilibrium ||, there is no rigorous way to determine whether the 
electroweak symmetry is restored non-thermally or not. Nevertheless, we know that the 
damping term causes the amplitude of the inhomogeneous modes to decrease according to 
Eq. (pO"D, and thus if there is no Higgs condensate present we expect |0| 2 to decrease as 
exp(— yt). This is indeed what we observe until t w 0.8 GeV -1 . We can therefore conclude 
that the electroweak symmetry is effectively restored during this time. 

At t>0.8 GeV -1 , (\(p\ 2 ) starts to grow and approaches its vacuum expectation value. 
Although it is not entirely clear that our approximations are still valid at this time, because 
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FIG. 3. The evolution of the renormalised \<f)\ 2 defined in Eq. (19) as a function of time. Because 
of the inhomogeneous fluctuations, \(f>\ 2 never vanishes, but its exponential decay shows that the 
electroweak symmetry is effectively restored up until i>0.8 GeV -1 



£>7 -1 , conservation of energy implies that the system must eventually end up in the bro- 
ken phase, because the initial energy density is lower than that at the transition point in 
equilibrium. 

In order to obtain a better insight to the dynamics of the fields we also measured the 
power spectrum of the electric field E{ of the U(l) gauge field B^. It is a gauge-invariant 
quantity, and although Gauss's law Eq. ([[(]) fixes the value of its longitudinal component, 
it does not affect its transverse components. In thermal equilibrium at temperature T, each 
Fourier mode Ej-(k) of the transverse components satisfies 



d 3 k 
(2tt)* 



\Ej(k)\ 2 ) = 2T, 



(25) 



and therefore one can define even in a non-equilibrium setting an effective temperature for 
each Fourier mode separately [15| 



(26) 



In Fig. ||, we show the effective temperature of different Fourier modes of Ei at the 
beginning of the simulation and at time t = 0.1 GeV -1 . Note that although the electric field 
is initially in vacuum, its effective temperature is nonzero, T e g(k) ~ k/2. This is caused 
by the initial "quantum" fluctuations (|TT|), and corresponds to the vacuum energy lj/2 of a 
quantum harmonic oscillator. In fact, in a quantum system, the effective temperature would 
be related to the occupation number as 



n{k) 



Teff(fc) 



1 

2' 



(27) 
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At t = 0.12 GeV -1 , the long-wavelength modes with k<k*~ 200 GeV have reached a 
high temperature T e g m 10 4 GeV, but the short-wavelength modes are still in vacuum. (The 
evolution of the power spectrum at large k is due to the damping factor, but these modes 
remain vacuum modes as the counterterm is also decaying.) 

This form of the power spectrum justifies the classical approximation, because the oc- 
cupation number of the long-wavelength modes is high, n ~ 100, and they are the relevant 
modes for the time evolution. The main effect of the short-wavelength modes is to renor- 
malize the couplings, and that has been taken into account perturbatively (|16|). 

As argued in Ref. the cutoff momentum is parametrically of order fc* ~ g<f) Q , and 
the effective temperature T e g ~ <fio/9- This is compatible with the observation |^9| that the 
theory becomes non-perturbative at momenta k < g 2 T, and thus the interactions below 
are not suppressed by the coupling constants and thermalization is more effective. 



Eventually, k* moves towards higher momenta, but this process is relatively slow |L5 
Instead, the thermalization is dominated by the decays of the long-wavelength bosons into 
fermions, which is approximated by the damping term. 

Obviously, it would be more interesting to follow the evolution of the power spectrum 
of the SU(2) field, but because its electric field is not gauge invariant, this is not possible. 
There are indications JJIJ that in the "vacuum" configuration flTSP, the cascade of power 
from the ultraviolet to the infrared modes is significantly faster in the SU(2) field than in 
scalar or U(l) fields. This would lower the ultraviolet part and raise the infrared part of 
power spectrum in Fig. |4j. However, since the power spectrum is already dominated by 
the infrared modes, this extra contribution to the infrared is not likely to be significant. 
Likewise, once the infrared modes have been excited, their occupation numbers are so high 
that the precise form of the ultraviolet power spectrum is unimportant. 

It is also possible that the energy transfer to the opposite direction, from infrared to 
ultraviolet, is faster in the SU(2) case, which would mean that the whole SU(2) field could 
thermalize much faster than the U(l) field in Fig. Even if this is the case, it would not 
presumably change the qualitative behaviour, because the SU(2) gauge field only constitutes 
a small fraction of the total field content of the standard model; full thermal equilibrium is 
not reached before all the fermions have equilibrated as well. Therefore we believe that our 
classical approximation gives the correct qualitative picture of the dynamics. 

During the period of symmetry restoration, baryon number is not conserved, and since 
the effective temperature of the long-wavelength modes is T e g- ~ 10 3 GeV, we can estimate 
using Eq. (ID) that the sphaleron rate is roughly 



sph 



10 6 GeV 4 . (28) 



Thus the non-equilibrium oscillations ought to be able to generate baryon asymmetry ||. 
On the other hand, the rate at which the temperature decreases, 7 ~ 2 GeV, is much higher 
than the baryon decay rate 7^ ~ 0.01 . . . 0.1 GeV in Eq. (|5]), which means that the sphaleron 
rate drops to zero so rapidly that the baryon asymmetry should not be washed out. 

VIII. EVOLUTION OF THE CHERN-SIMONS NUMBER 

In thermal equilibrium, the sphaleron rate r sp h can be measured directly, because Nqs 
performs a random walk and r sp h is given by the diffusion constant 
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FIG. 4. Effective temperature (26) of different Fourier modes of the hypercharge field measured 
in the initial configuration and at t = 0.1 GeV -1 . At this time, the Higgs fluctuations have heated 
up the long-wavelength modes to a high effective temperature. For comparison, the plot also shows 
the effect of the damping term to the vacuum, measured in a run with 0rj = 0. 



sph 



lim 

t— >oo 



((N cs (t) - N cs (0)) 2 ) 



t 



(29) 



In Ref. m, this was generalized to non-equilibrium cases by defining the time- dependent 
sphaleron rate 



T sph (t) = j t ((N cs (t)-N cs (0)f). 



(30) 



This was then used in a phenomenological Boltzmann equation to estimate the baryon 
asymmetry generated during preheating. 

In Fig. |5|, we show the time evolution of the variance, ((iVcs(i) — -/Vcs(O)) 2 ), of the Chern- 
Simons number during the simulation (upper plot). Instead of growing monotonically, as 
Eq. (|30| ) assumes, the variance oscillates until t ~ 0.15 GeV -1 , and then finally settles down 
to a value that is well below its maximum during the oscillations. This indicates that the 
behaviour of Nqs is not purely Brownian, but it contains a deterministic component, which 
tries to decrease its value. To see if this is a result of the damping we performed simulations 
where 7 = 0. The lower graph in Fig. |5] shows that the oscillations persist even if we remove 
the damping term. In that case, Nqs keeps on wandering after the initial oscillations because 
the system thermalizes to a high temperature, but the behaviour is still non-Brownian at 
t ~ 0.1 GeV -1 . In fact we would not expect the damping to have an effect as the damping 
timescale is 0.5 GeV -1 and the effect occurs before 0.1 GeV -1 . 

This deterministic behaviour is explained by the interaction between the gauge and 
Higgs fields. When the Higgs field <£> is non-vanishing, one can define the integer-valued 
Higgs winding number Nh [|l|] • 
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FIG. 5. (Nq S ) for dissipation 7 = 2 GeV" 1 (top) and 7 = (bottom) with a cooling depth 
of 0.9a 2 , averaged over 12 runs. Initially (Nq S ) oscillates and then starts to decrease vanishing 
eventually, which means that no topological transitions actually took place. In thermal equilibrium, 
(Nqo) would grow linearly, which is what happens in the dissipationless run at late times. 



N. 



H 



d 3 xe ijk Tr (di&djQdk&Q 



24vr 2 J " " V % * V • ^ 31 ) 

In particular, in a vacuum iVcs = Nh, because any vacuum configuration can be gauge 
transformed into one in which the gauge field vanishes, which implies iVcs = 0, and the 
Higgs field is constant, which implies Njj = 0. 

In previous studies of iVcs i n thermal equilibrium |32| , |33| , |3^ , |35| , |36| , |37|| , this has not been 
a problem. If we write 

N C s + N H Ar N cs - N H 



N, 



cs 



N + + N_, where A^ 4 



(32) 



2 2 

the different vacua are labelled by N + . In thermal equilibrium, N + performs a random 
walk because of sphaleron transitions, and iV_ merely fluctuates around zero. In our non- 
equilibrium case, these fluctuations must die away as the system cools and therefore iV_ 
must tend to zero. Thus, iV + will be the sole contributor to iVcs in the final state. Given 
that, the contribution to (Nq S ) from iV + will be monotonic and all the early oscillations of 
Fig. |5| are due to iV_. This implies that practically no topological sphaleron transitions took 
place in any of our runs, while Eq. ( p8|) would have predicted around 5000 of them. 

We can further back up this claim that the changes in (Nq S ) are not of a topological 
nature by looking at (n^g), the U(l) equivalent of (Nq S ), [see Eq. @]; large gauge trans- 
formations leave (ncs) unaffected and so its evolution is due only to fluctuations. If there 
is no topology change in the SU(2) sector then we can relate the fluctuations of (A^q S ) to 
(7i CS ) via II 



cs)/<"cs> - Hg/g'f 



(33) 
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FIG. 6. The squared U(l) Chern-Simons number (n^ s ) multiplied by the factor [3(g/g') 2 ] 2 
[see Eq. (|33| )]. It is of similar order of magnitude as (Nq S ) in Fig. [|, which supports the idea that 
(Nq S ) consists of fluctuations only. 

where the factor 3 is because there are three SU(2) vector bosons to the one U(l) vector. 
The plot of [3(g / g') 2 ] 2 (uqs) i n Fig- || shows that the U(l) Chern-Simons number is of the 
order expected if the evolution does not contain topological changes. 

This non-Brownian element in the motion of Nqs means that the simple Boltzmann 
equation does not give an adequate description of the dynamics. The same effect is present 
also in the CP violating case, and this may lead to baryon washout even well below T c if the 
Higgs winding remains small during the non-equilibrium time evolution. Even if the initial 
fluctuations generate large Nqs, and therefore a large baryon asymmetry, the interaction 
with the Higgs field may destroy it soon afterwards. 

IX. CONCLUSIONS 

In this paper, we have studied the out-of-equilibrium dynamics of the electroweak theory 
in a case where most of the energy is concentrated in the long-wavelength modes. These 
initial conditions are interesting because they resemble the state after preheating at the elec- 
troweak scale, where the Higgs has a large initial value. All the other degrees of freedom were 
in vacuum, and were given initial conditions that reproduce the tree-level two-point function 
of the quantum vacuum. The coupling to the fermions was approximated by a phenomeno- 
logical damping term, whose effects were also taken into account in the renormalization by 
using a time-dependent mass counterterm. 

Our method can be used to simulate more realistic scenarios of electroweak preheating 
in a straightforward way by coupling the Higgs to the inflaton field. In order to achieve 
accurate results, one would also have to be able to improve on our approximations, but 
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we believe that even in its present form, the method can give reliable order-of-magnitude 
estimates. 

Our results show that the oscillations of the Higgs field can lead to non-thermal restora- 
tion of the SU(2) gauge symmetry. This supports the scenario of baryogenesis where the out 
of equilibrium state, required by the Sakharov criteria, is provided by a form of parametric 
resonance rather than a thermal phase transition but a truly quantitative test would 
require a measurement of the generated baryon asymmetry. In our simulations this was not 
possible because the equations of motion were CP invariant. 

However, we also studied the Chern-Simons number of the SU(2) gauge field and showed 
that its motion is not pure random walk. We believe this is due to interactions with the 
Higgs fields which force the vacuum value of the Chern-Simons number to equal the winding 
number of the Higgs field. If the Higgs winding therefore does not change significantly, 
these interactions will reduce Acs, destroying the baryons previously created. Conversely, if 
a mechanism greatly enhanced Nh then baryons would be created as the winding numbers 
equilibrated. In every one of our simulations, the final Chern-Simons number vanished, which 
suggests that the rate of true topological transitions that would lead to stable baryons is 
lower than expected during the non-thermal symmetry restoration. It remains to be seen 
if this is true in more realistic models as well. In any event, it means that the baryon 
asymmetry cannot be estimated using any effective description that does not take this effect 
into account. 
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APPENDIX A: EVOLUTION EQUATIONS 



To be explicit we give here the lattice method p2| , p5| used to simulate the continuum 



gauge theory. The starting point is the gauge covariant derivative on a lattice, 

D4 = - (Ui(t, x)Vi(t, x)<f>(t,x+i) - (pit, x)) , (Al) 
a 

D <P = - (U (t, x)V {t, x)<j>{t + a u x) - (f>{t, x)) . (A2) 
at 

Here the ~ refers to the lattice covariant derivative, and a and at are the lattice spacings in 
the space and time directions, respectively, with x+i referring to the position one lattice site 
away in the i direction. The fields Ui(t,x) and Vi(t,x) are fundamental SU(2), U(l) valued 
respectively, and correspond to the continuum gauge fields according to, 

U i (t,x)=exp(~gaWta a ), (A3) 
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U (t,x) = exp ( --ga t W£a L 



Vi(t,x) = exp (^—g'aB^j , 
V (t,x) = exp (-^g'atB^j . 
One then finds that the standard continuum gauge action is modelled by 



s = e a t aZ 

x,t 



i 

2 ETr(^)--^ETr(^) 



(A4) 
(A5) 
(A6) 

(A7) 



a 2 t a 2 g 2 



a 2 a 2 g' 2 



Here we do not use the repeated index summation convention. The field strengths for the 
gauge fields are given by 



Uij(t,x) = Uj^x^^x+Dtf^x+^Uj^x), 
Vij(t,x) = VjfaxWfax+fiVjfax+ilVlfax), 

which means that the lattice action has the following gauge symmetry 

4>(t,x) — > Qi(t,x)Q2(t,x)(f)(t,x), 
Ui(t,x) — ► D,2(t,x)Ui(t,x)£V 2 (t,x+i), 
Vi(t,x) — ► fii(t,a;)\/ i (t,x)fi t 1 (t,x+?), 
Vl 2 G SU(2), G U(l). 



(A8) 
(A9) 



(A10) 
(All) 
(A12) 
(A13) 



In our simulations we use the equivalent of temporal gauge, in which the link variables Uq 
and Vq assume the identity value, with their equations of motion leading to Gauss type 
constraints. The equations of motion are found by the lattice equivalent of functional dif- 
ferentiation of the lattice action, 



4>(t+a t , x) = (p(t, x) + a t ir(t+a t /2, x), 
Vi(t+a t ,x) = ^g'aa t Ei(t+a t /2,x)Vi(t,x), 

Ui(t+a t ,x) = gaa t Fi(t + a t /2,x)Ui(t,x), 
7i(t+a t /2,x) = (l-aa)n(t-a t /2,x) 



(A14) 
(A15) 

(A16) 
(A17) 



+a t 



- E (Ufa x)Vi(t, x)<j>(t, x+i) - 2<f>(t, x) 



lm(E k (t+a t /2,x)) 



+Ul(t,x)V?(t,x)<P(t,x-i)) - 
[l-a tl )lm{E k {t-a t /2,x)) 



dV 

90t 



(A18) 
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+a t 



lm[tf(t, x+k)Ul(t, x)V k \t, x)(j)(t, x)} 



—r-n J2 lm l V k(t, x)Vi(t, x+k)V k \t, x+i)V?(t, 
g'a A Y 



x 



+Vi(t, x-i)V k {t, x)V?(t, x+k-i)V k \t, x-i)] 
Tr(i<j m F k (t+at/2,x)) = (l-a a )Tr(ia m F k (t-a t /2,x)) 



(A19) 



+a t 



ga 



Re[(j)\t, x+k)Ul(t, x)V k \t, x)ia m (j)(t, x)] 
^ E Tr [ i(jmU k (t, x) Ui (t, x+k)Ul(t,x+ 1) U} (t, x) 



+ia m U k (t, x)U}(t, x+k-i)U\{t, x-i)Ui(t, x-i)} 



Here we have included a damping factor 7 which is a phenomenological term put in by hand 
to model the decay into fermions. Such a term does not spoil the gauge invariance of the 
model. The remaining components of E k and F k are found by using 



det(E) 



det(F) 



1 



g aa t gaat 
The Gauss constraints which are implied by the lattice gauge symmetry are 

- hn[Ei(t+at/2, x) - E i (t+a t /2,x-i)] 
a i 

-g'lm[n ] (t+a t /2, x)<p{t, x)\ = 0, 

- Tr[iv k Ei{t + a t /2, x) - ia k Uj{t,x-i)F i (t+a t /2, x-i)Ui(t, x-i) 
a i 

-gRe[n\t+a t /2, x)ia 4>(t, x)} = 0. 



(A20) 



(A21) 



(A22) 



APPENDIX B: RESTORING THE GAUSS CONSTRAINTS. 



When the random initial conditions on the gauge and scalar fields are imposed we find 
that the Gauss constraints will generically be violated. In order that the simulation is to 
make sense this must be remedied. In the Abelian case the form of Gauss's law is simple 
in k space so the longitudinal component of Ej may easily be set to the relevant charge 
density. In the non Abelian case however the derivative in Gauss's law requires a gauge field 
to make it gauge covariant, it is not such a simple matter therefore to find the equivalent 
of the transverse and longitudinal components. The way we get around this is to relax 
the unphysical gauge degrees of freedom in the electric fields and scalar momenta into a 
state which satisfies Gauss's law. This is done by making a Hamiltonian out of the Gauss 



constraint and using this to define a dissipative motion |32| . [23 



Defining G° to be the hypercharge Gauss constraint expression ( A21 ) and G k to be the 
remaining three expressions for the SU(2) field ( |A22p , such that we require G° = 0, G k = 0. 
Now introduce the Hamiltonian 
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H = J2a 3 [G°G° + G k G k 



(Bl) 



By writing F { = Ff — ^a a F^ we then evolve the 7r, E h Fi fields according to the dissipative 
equations 



<9r ' 



where represent the generic momentum field, to find 

■ tg'G (T,x) ( f>(T,x)-tgJ2v k G k (T,x) ( j ) (T,x) 

k 

^) = -%V,*)-G W +j) 



dn 



dr 



dr 



&(t,x) - -Tr (a k Uj(r,x)a%(r,x)) G k (r,x+j) 



(B2) 



(B3) 
(B4) 
(B5) 
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